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Abstract. An ensemble of directed macromolecules on a lattice is considered, where 
the constituting molecules are chosen randomly with A'^ different colors. Molecules of 
the same color experience a hard-core (exclusion) interaction. We study the robustness 
of the macromolecules with respect to breaking and changing the color of a constituting 
molecule, using a 1/N expansion. The properties depend strongly on the density of 
macromolecules. In particular, the macromolecules are robust against breaking and 
changing the color at high densities but fragile at intermediate densities. 
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Figure 1. Two directed macromolecules. 

1. Introduction 

In a liquid of molecules long chain-like macromolecules (MM) can be formed due to a special 
steric intermolecular interaction. These MM are often called polymers 1 . We assume that 
the MM are directed, i.e. they have a preferred direction in the d-dimensional space. This 
implies that they are stretched and do not form loops or overhangs. 

In many systems the constituting molecules (CM) of the MM are identical. However, 
in more complex systems, like in biological substances, the CM are chemically different. An 
example for the latter are the sequences of amino acids in proteins O El El and molecular 
sequences along a DNA helix. For our statistical model we assume that there are N different 
kinds of molecules, represented by N colors. They are randomly chosen to contribute to 
the MM such that the resulting MM can be viewed as a chain of colored beads (cf. Fig. 
1). An ensemble of MM can be considered as being created in a liquid of molecules by a 
steric interaction among the CM. As a result, we have a space-filling system of molecules, 
where some are constituting the MM and others exist as individual molecules. Properties 
of the ensemble on large scales should not depend on local structures. This allows us to 
consider a lattice model 01, where each lattice site is occupied by a molecule of each type. 
It is assumed that molecules of the same type have a strong repulsive interaction. This 
restricts us to at most one molecule of each type at any given lattice site. The fact that the 
interaction of different molecule types is neglected can be understood as a chemical property 
of the molecules. 

Any pair of neighboring molecules can form a bond state, representing a pair of 
molecules, which acts as a building block for an MM. There is a tendency of these bond states 
to bond with other molecules. This leads eventually to an equilibrium state, where we have 
a mixture of only large MM and individual molecules. The repulsive interaction between 
the molecules creates correlations among the MM which lead to interesting properties of the 
ensemble. For instance, the density of MM as a function of the probability of a local bond 
state is of interest, as well as the question how likely it is to break an MM or to replace 
a CM of special type by one of a different type. The latter can happen during a collision 
process when an individual molecule collides with an MM and exchanges with a CM of the 
MM. 

We consider an equilibrium ensemble of lattice MM and define three quantities to describe 
its global properties with respect to the molecular structure. There are two quantities that 
measure the correlations in real space. One measures how likely it is to break a directed 
MM and take it apart in space (cf. Fig. 2), the other one measures how likely it is to 
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Figure 3. Changing tlie type of a constituting molecule inside a macromolecule, caused 
by a collision with a molecule from the surrounding liquid. 



create a new (finite) MM by clianging the type of the CM in an MM. If the ends of the 
broken MM or the ends of the new (finite) MM are x and x' we write for the corresponding 
quantities Gxx' and Hxx', respectively. It will be shown that Gxx' and Hxx' are directly 
related. Another quantity is the coherence in a (molecular-component) space. This is called 
Saa' and measures how likely it is to replace a CM of type a in an MM by one of type a' 
(cf. Fig. 3). All these quantities are evaluated in the limit N ^ oo and in a expansion. 
We apply a functional-integral formalism developed in Ref. |S] to perform the calculations 
in a straightforward manner. 

2. The Model 

A statistical model for the equilibrium statistics of directed MM was introduced in Ref. 

We briefly summarize the construction of the model and then turn to the evaluation 
of the interesting correlations mentioned in the Introduction. The model is based on the 
idea that on each site of a (cubic) lattice N different molecules are located. They exist 
either as individual (separated) molecules of a liquid with relative weight fi or contribute to 
a directed MM as a bond state of a nearest-neighbor molecular pair at x, x' of type a and 
a'. Moreover, the MM are directed, i.e. there is a special direction t in the d-dimensional 
space as a preferred direction of the MM. The fluctuations of the MM are perpendicular to 
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this, indicated by the (d — l)-diniensional vector r. Therefore, we write x = {t,r) for a site 
on the lattice. The relative weight of a nearest-neighbor molecular pair at x, x' of type a 
and a' is 



1 \ J/N a t' — t + 1, r — r' or r, r' nearest neighbors 



N ' I otherwise 

The directed chain-like structure is given by the matrix w that is asymmetric in the direction 
t but symmetric with respect to r. x' is assumed to be 'later' along the direction of the MM. 
w is also symmetric with repect to a,a' . Thus a and r are of the same type in contrast to 
t. 

The statistical distribution of configurations of MM has to take into account a strong 
repulsive interaction between molecules of the same type at the same site. This can 
be conveniently be described by a commutative algebra of nilpotent variables 77"'°^ (i-e., 
{Vx''^y = if Z > 1) with a — 1,2,.. .,N and a — 1,2 5 . Each molecule, characterized 
by its position x = {r,t) on the lattice and by its chemical character (a, a), is described 
by a variable. Since {jy"'"^} are nilpotent, the repulsive interaction results in an exclusion 
principle for molecules of the same type. The atoms are considered as classical objects. 
Consequently, products of rj"''^ must be commutative. A segment of the MM (i.e. a bond 
state of two molecules or a dimer) reads 

ckck' ck 1 ?2 
'^x.x'Vx' Vx' ■ 

A linear mapping of the algebra {tj"''^} to the complex numbers can be written as an integral 
on a lattice A with 6 



n n n 



if L'^L, I, = 1,2,...,7V, j,,„ = l,2 
otherwise 



The integral vanishes if the product is incomplete with respect to the lattice A, to 
a — 1,2,.. .,N or to CT = 1,2. It can be used to write the statistical weight Pj of a 
configuration / of MM as 

with 

(x,Q;x',a')G/ (x^a) 

where (x, a; x' , a') is a bond between two CM of type a and a' at sites x and x' . The 
normalization Z is the sum over all possible configurations {/}: 



{1} 

Finally, we impose periodic boundary conditions on the lattice. This implies a system on a 
torus with closed loops of MM. However, the choice of boundary conditions should not be 
crucial for the local properties of the MM like the density of the MM, Gxx', Hxx' ^ or Saa' ■ 

3. Functional-Integral Representation 

Introducing analytic functions of the nilpotent field and using the properties of the integral, 
we may write for the normalization 

f f ^ 1 

^ = E 1^1= /<5xp[E E ij^w^,^' + l^^x.x'Sa,a')r]^-^r]'^,'^ . (1) 
{/} x^x' Q,a' — 1 
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N ^ 

{/} x,x' a,a' — 1 

in Z . Another important consequence of the properties of the nilpotent field is that we can 
calculate the probability that a point a;, a belongs to a MM configuration /: 

If, , { Pi if (x, a) e / 



Z J XX 1 Q otherwise 

After summation over all configurations {/}, the probability of a molecule of type a at site 
X being part of any MM is 

{1} {1} 

This is also the local density of an MM. In a similar manner, we can express the 

correlations in terms of normalized integrals 

(...) 4 /...»' 

as correlations in real space 

Gxx' = {Vx"^Vx'^) = -C{xa;x'a), 
Hxx' = iVx^Vx'') = -C{x'a;xa), 

and correlations in a space 

5aa' = (77"'S"''') = -C(xa;a;a'). 

These are correlation functions in real space and in a space, respectively. The density of 
MM can also be expressed as 

n(/i) = 1 — ^C{xa; xa). 

A lengthy but straightforward calculation 5 shows that the normalization Z in Eq. 
can be written as an functional integral with respect to complex fields and x 3-8 

Z= j e-'^'^Wd^xdXx. (2) 

•' X 

with 

5" = (1 + + (X, X) - X! + ("^^ + iXx){(px + iXx)]- 

X 

An analogous calculation for the correlation functions leads to 

^■■■) = i /■■■ e-''^\{d^xdXx. 

X 

The correlation functions then become 

C{xa; xa) = (\^n+{ipx+ iXx){,V*x + %)]~^) 

for (x, a) ^ (x', (3). 



[^i + {ipx + iXx){^*x + iX*x)][l^ + {^x' + iXx'){^Pl, + ix*x-)] 
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3.1. 1/N Expansion 

The integrals in Eqs. Q and © depend on the number of different molecule types N 
only through the prefactor N in the exponential. This enables us to perform a saddle-point 
integration for large values of iV [S]. The result is a 1/N expansion. The homogoneous 
saddle point of the integral in Eq. ^ satisfies the relation 

1 _ 1 _ 

Then there is a trivial solution (ySo = <^o = Xo = Xo = and a non-trivial solution with 

ipii^i = 4(1 - ^). 

An expansion of the exponent in Z in terms of the deviations around the saddle point 

(t>x = + ix + ^(f'x with 54>3: -.^ Sifix + iSxx, ^'i>*x ^^*x + "^^xl 

provides corrections in powers of 1/N . For a non-zero density of MM (i.e. C = /^^^ > l), an 
expansion up to second order leads to 



+Va^mx) + (c - If mi?) + (c - i?{hihx') - (c - mKHi') 

'{C - i)^aC-i){Hx) + (C - ifiHl) - 2(C - 1) - (C - i)VC(C-i)W.') 



-{c-i)mi,f) + ./a^mi,)], 

The first term is the large- limit, all other terms either vanish because of {5(j)x) — or 
are of order because of (Scj)-^,) — o{l/N). The saddle-point calculation leads to the 
result that the first order terms vanish, and the second order terms can be expressed by the 
integrals 

-r',t- t') {Scf^Jc^x,) = {S^Ml') « 



■ e 



N J Jo 2B(fc)cos(w) -C + (C-2)B(A;)2 27r (27r)''-i 

$(r-r',t-t') {S^Ml') ~ 
1 f f^^ B{k){B{k) -Ccosjcj)) d'^-'k 



N J Jo 2B(fc)cos(tj)-C + (C-2)B(fc)2 2tt {2ttY- 

with the function B{k) := 1 — J + -j^ cos kj. 

4. Discussion of the Results 

The density of MM is 



7l(/i) 



1 - ^ + o{l/N) for < 1 
for II > 1 



It increases with decreasing N, as shown in Fig. 4. This means that it is more likely to form 
a MM if the number of available types of molecules is reduced. 

How likely it is to break a MM depends on the position of the remaining ends x — x'. The 
corresponding quantity Gxx' decays with |a; — a;'| which indicates that the molecules have a 



Correlations in Systems of Complex Directed Macromolecules 



7 




Figure 4. The density of macromolecules and the correlations for breaking a 
macromolecule (G) and replacing a molecule in a macromolecule (5). 



tendency to form a MM. However, it does not vanish asymptotically but stays nonzero for 
C > 1 and = t: 

G= lim C(xa;a;'a)«(l-C"^)(C"^ + 

The C dependence of this asymptotic result is shown in Fig. 4. The same property is valid 
for the creation of a finite chain of molecules with end points x,x' , decribcd by Hxx' ■ 

The correlation in a space Saa', describing how likely it is to change the type of 
molecules in the MM, does not depend on the specific pair a, a' as long as a 7^ a' because 
the model connects all a ^ a' with equal probability: 

S = C{xa; xa') « ("^ - + J (l - ^ + ^) *(0, 0) + ^ (l - ^) '^(0, 0) 

The N ^ 00 terms of both correlations G and S agree but there is a difference in terms of 
order 1/iV (s. Fig. 4). It is unlikely to break an MM and separate the pieces at low and 
at high density. This is a consequence of the fact that there are either too few MM or too 
many at high densities such that the separation of the pieces is difficult. It is remarkable 
that already at densities less than 1 the separation is completely blocked, as indicated by 



U(o,o) ). 
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the vanishing G ior < n < 1. Moreover, the substitution of a CM is most likely at 
intermediate densities but suppressed at low and high densities. Again, there is a critical 
density < 1 such that the substitution is blocked for < n < 1. This can be understood 
as an effect of interaction: Since there is a exclusion between molecules of the same type, we 
can only substitute if there is no contribution of this molecule from another MM. At high 
densities it is very unlikely that a site is not occupied yet by a CM of a specific type from 
another MM. 



4.I. Algebraic decay of the correlation function 



To analyse asymptotic behaviour of the correlation function for large distances 



of a 



three-dimensional system, we first have to perform the w-integration in Eqs. H3.1|) and H3.1|l . 
This leads to 



1 

$(r-r') = — 



-B(fc)(C-l) 



d^k 



(27r)2 ' 



l>(r - r') 



1 

TV 



C 



2 V((C-2)S(fc)2-CP-4i?(fc)2 
The spacially decaying part of the correlation function is given as 



(27r)2 



2 



1 



For small values of k we can assume B{k) w 1 in the numerators of <l>(r — r') and $(r — r') 
The square root in the denominators can be approximated by using B{k) = l-^ + 0(fc4). 



V((C-2)i?(fc)2-C)2-4B(fc)2 = V4(C - + 0(fc4) « 2^/C^\k\ 
Therefore, the asymptotic behavior of the correlation function is proportional to 



const 



(2^)2 



OC const 



27r noo 



\r — r'\ ^ const 



"'0 



where k has been substituted by k — fc|r — r' 
l/|r-r'| ford = 3. 



=const 

Thus the spatial correlations decay like 



5. Conclusions 



A statistical ensemble of directed MM, constructed with a random sequence of N different 
molecules, has been treated for large values of N in terms of a 1/N expansion. We have 
studied the robustness of the MM with respect to breaking and replacement of individual 
CM. Breaking a MM and taking the pieces apart in space reveals a long-range correlation 
with an algebraic decay with the inverse distance. The properties in space and in a are 
identical in the N ^ oo limit but different in 

At high densities the broken MM cannot be separated to arbitrarily large distances due 
to the interaction between different MM. The replacement of individual CM is most likely at 
intermediate denstities but strongly suppressed at high densities. This is also a consequence 
of the repulsive interaction between the molecules of the same type. 
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